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PART  II:  THE  ZERO  PRESSURE  GRADIENT 
LAMINAR  BOUNDARY  LAYER 


by 

John  0.  Powers  and  Edgar  Krahn 


ABSTRACT:  It  has  been  previously  shown  (reference  (a)),  that 

the  two-thickness  integral  method  of  determining  boundary- 
layer  characteristics  provides  reasonable  results  for 
stagnation-point  flows  even  when  the  effects  of  equilibrium 
air  dissociation  are  considered.  This  result  implied  that  the 
method  may  have  more  general  application  and  it  therefore 
seemed  advisable  to  explore  additional  characteristics  of  the 
method  in  order  to  further  assess  the  extent  of  its  applica¬ 
bility.  The  zero  pressure  gradient  case,  with  its  associated 
limiting  cases,  i.e.,  zero  Mach  number  and  zero  heat  transfer, 
has  been  used  as  a  basis  for  this  further  evaluation  primarily 
because  of  the  availability  for  comparison  of  computations 
performed  by  "exact"  methods.  In  general,  it  appears  that  of 
the  profile  representations  considered,  a  two-parameter,  sixth- 
degree  velocity  profile  gives  the  most  desirable  results  for 
a  flat  plate  when  used  with  the  present  method. 
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local  skin-friction  coefficient 
mean  skin-friction  coefficient 
total  enthalpy  thickness 
velocity  ratio,  u/U 
total  enthalpy  ratio,  H/H^  -  1 
static  enthalpy 
recovery  enthalpy 
total  enthalpy 

polynomials  used  in  equation  (24)  where 
i  “  1,  .  .  .  .9  (coefficients  of  polynomials 
given  in  Table  I) 
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static  pressure 
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heat-transfer  rate 
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temperature 

X  direction  velocity  component 
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y  direction  velocity  component 

R  direction  velocity  component  defined  by 
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co-ordinate  along  body  in  stream  direction 
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and  (25)  where  1  “  0,  1,  8  (given  in  Table  I) 

04  coefficients  given  by  equation  (36)  where 

i  =  1,2 


6 

6* 
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viscous  boundary-layer  thickness  in  x.q  plane 

displacement  thickness  in  x,q  plane 

thermal  boundary-layer  thickness  in  x,q  plane 
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first  velocity  profile  form  factor  defined 
by  equation  (11) 

second  velocity  profile  form  factor  defined 
by  equation  (11) 

enthalpy  profile  form  factors  defined 
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maximum  ((5,A) 
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HEAT  TRANSFER  IN  DISSOCIATED  AIR  BY  A 
TWO-THICKNESS  INTEGRAL  METHOD 
PART  II:  THE  ZERO  PRESSURE  GRADIENT 
LAMINAR  BOUNDARY  LAYER 


INTRODUCTION 

1.  The  two-thickness  integral  method  of  computing  laminar 
boundary-layer  characteristics  at  the  stagnation  point  of 
bodies  including  the  effects  of  dissociated  air  in  eouilibriura 
has  been  presented  in  reference  (a).  From  that  investigation, 
it  appeared  that  the  method  may  have  more  general  application 
and,  accordingly,  it  was  considered  desirable  to  investigate 
additional  characteristics  of  the  method.  The  procedure  pres¬ 
ented  was  based  on  an  extension  of  the  Karman-Pohlhausen  method 
as  reported  by  Morris  and  Smith  (reference  (b)),  and  had  been 
modified  to  include  the  latest  information  on  the  thermodynamic 
and  transport  properties  of  air  as  determined  by  members  of  the 
Thermodynamic  Section,  National  Bureau  of  Standards  (NBS) .  The 
further  investigation  has  been  conducted  on  the  zero-pressure 
gradient  boundary  layer  and  the  results  havv?  been  compared  with 
the  results  of  "exact”  methods  (references  (c),  (d) ,  (e),  and 
(f)). 

2.  Van  Driest  (reference  (d))  and  Klunker  and  McLean  (reference 
(e)),  both  consider  the  temperature  variation  of  the  fluid 
properties  but  do  not  perform  calculations  for  temperatures  at 
which  air  becomes  dissociated.  The  property  values  used  in 
references  (d)  and  (e)  were  different,  hence,  discrepancies' 
between  these  results  may  be  attributed  primarily  to  this 
difference.  The  method  of  solution  of  the  basic  differential 
equations  was  similar  in  that  a  transformation  was  used  to 
obtain  an  equation  in  integral  form  which  in  turn  was  solved 

by  the  method  of  successive  approximations.  The  transformations 
were  different,  however,  because  Van  Driest  used  Crocco's  pro¬ 
cedure  which  results  in  velocity  as  the  independent  variable 
whereas  Klunker  and  McLean  used  the  procedure  of  Schuh  which 
yields  normal  distance  as  the  independent  variable. 

3.  The  recent  work  of  Wilson  (reference  (f))  which  is  based 
on  the  Crocco  procedure,  perm.its  comparison  with  the  two- 
thickness  Integral  method  in  temperature  regions  where  dis¬ 
sociation  affects  the  property  values.  Wilson  used  the  most 
current  NBS  thermodynamic  and  transport  data  available  at 
the  time  of  his  investigation;  however,  since  the  amount  of 
transport  property  data  was  limited  in  the  dissociation  region, 
it  was  necessary  for  him  to  devise  an  interpolation  scheme 
which  produced  reliable  results  over  a  wide  range  of  stream- 
to-reference-density  ratios.  Since  the  completion  of  the 
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reference  (f)  investigation,  the  NBS  (reference  (g))  has 
generated  new  transport  property  data  in  the  temperature  range 
between  2700OR  ;and  6000°^*  The  two  sets  of  property  data 
compare  very  favorably  with  the  largest  discrepancy  (approxi¬ 
mately  five  percent),  appearing  in  the  Praadtl  number  at  280O‘^R. 
The  reference  (g)  data  have  been  incorporated  herein  and  the 
interpolation  scheme  of  reference  (f)  has  been  used  for  viscosity 
determination  at  temperatures  above  6000OR. 

4.  In  Part  I  of  this  investigation  (i.e.,  "Stagnation  Point 
Heat  Transfer,”  reference  (a)),  a  one  parameter,  fourth-degree 
polynomial  velocity-profile  representation  was  used  because 
of  its  known  suitability  for  stagnation  point  flows.  When 
applied  to  general  cases,  this  representation  results  in  a 
lisnitation  because  of  the  single  parameter,  namely,  that  its 
absolute  value  cannot  exceed  12,  Such  a  limitation  is  charac¬ 
teristic  of  methods  based  on  the  Pohlhausen  procedure.  Ac¬ 
cordingly,  in  the  Part  II  investigation  reported  herein  a 
two-parameter  velocity-profile  representation  was  used  in  an 
attempt  to  extend  the  range  of  applicability  of  the  method. 

In  addition,  foxirth-,  fifth-,  and  sixth-degree  velocity  poly¬ 
nomials  have  been  used. and  total  as  well  as  static  enthalpy 
profiles,  botli  of ‘the  sixth-degree,  were  investigated.  While 
altering  the  velocity  profile  representation  avoids  some 
limitations,  other  limitations  are  introduced.  The  zero 
pressure-gradient  case  has  the  advantage  of  presentin..  an 
opportunity  to  explore  some  of  these  limitations  with  at 
introducing  excessive  complexity. 

ANALYSIS 

5.  The  complete  equations  of  motion  including  pressure 
gradient  for  the  flow  of  a  real-gas  compressible  fluid 

given  In  reference  (a) .  For  a  two-dimensional,  zero  pressure- 
gradient  boundary  layer  the  equations  become: 


(y 


3i\ 


/  \  //  y 
cy^ 


r)  /  ^ 


O', 


u  ■ 


d 


(continuity)  (1) 


(momentum)  (2) 


(energy)  (3) 
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It  is  noted  that  the  Prandtl  number  used  in  the  energy 
equation  includes  contributions  of  both  the  conductive  and 
ordinary  diffusive  energy  fluxes  for  air  in  equilibrium  dis¬ 
sociation  . 


6.  As  in  the  reference  (a)  development,  the  partial  differ¬ 
ential  equations  (1),  (2),  and  (3)  are  initially  transformed 
from  the  x,  y  plane  to  the  x,  q  plane  by  the  Dorodnitsyn  trans¬ 
formation  in  the  form  y 

7:j  g--"/ 


As  a  resiilr  of  this  transformation,  the  continuity,  momentum, 
and  energy  equations  become: 


~f-  : -  -  Cd' 


w^£r/^E: 


Pi/ 

/  ■-_  u,  — ^  ~ 


(4) 


U  —  -L 


X  ^  =  - - ^  - 

ax  4?^  \'>5 


^  LL. 


(5) 


3>:.  '  c!?  "  e' 0?-^ 


(6) 


7.  When  it  is  assumed  that  the  total  enthalpy, 


/  / 


A 


;? 


at  the  edge  of  the  boundary  layer  is  constant  and  equations 
(5)  and  (6)  are  integrated  from  the  wall  to  the  outer  edge 
of  the  boundary  layer,  the  result  is: 


and 


c/e  Tr^ 

^  PU^ 

O 


PL 

LX 


(7) 


(8) 
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where  and  are  respectively  the  heat  transfer  and  shear 
stress  at  the  wail.  The  corresponding  integral  equations  for 
the  pressure  gradient  case  are  given  in  reference  (a).  The 
boundary- layer  momentum  thickness,  6,  and  total  enthalpy 
thickness,  E,  are  defined  in  the  transformed  plane  as: 


where  <5  and  ^  are  the  values  of  q  at  the  outer  edge  of  the 
viscous  and  thermal  boundary  layers,  respectively.  The  dis¬ 
placement  thickness,  <5*,  and  the  density  thickness,  ^ 
which  supply  additional  information  but  are  not  essential  for 
the  solution  of  the  system  of  equations,  are  defined  as: 


The  inverse  traiisf ormations  of  these 
the  viscous  and  thermal  thicknesses, 


e- 


thicknesses, 
are : 


[a^-J 


(9b) 


including 


10; 


where  quantities  in  the  physical  plane  are  indicated  by  the 
subscript  y  and  %  is  the  maxlraum  of  6  or  A, 


4 


NOLTR  61-25 


PROFILE  REPRESENTATION 

8,  The  foregoing  development  differs  from  the  Part  I  (refer¬ 
ence  (a))  development  only  in  the  fact  that  the  pressure 
gradient  terms  have  been  omitted  for  the  basic  equations.  In 
the  Part  II  Investigations,  an  attempt  was  made  to  extend  the 
range  of  applicability  of  the  method  and  accordingly,  a  two- 
parameter  velocity-profile  representation  was  used  instead  of 
the  one -parameter  representation  as  used  in  Part  I.  It  was 
also  considered  desirable  to  utilize  fourth-,  fifth-,  and 
sixth-degree  polynomials  for  the  velocity  profiles  to  evaluate 
the  relative  advantages  of  each  by  comparing  the  results  with 
other  methods.  Parametrically,  the  velocity  profiles  are 
presented  in  the  form: 


f 7)  =fJl^  *  \  <'7)  ’■  7-7) 


where  ti  “  n/d  and  the  parameters  and  X2  are  defined  by 
equation  (17)  .  In  forming  the  profiles,  it  was  considered 
desirable  to  utilize  as  many  wall  boundary  conditions  as 
practical  since  the  phenomena  of  interest  occur  primarily 
at  the  wall.  The  boundary  conditions  used  varied  with  the 
degree  of  the  polynomial  and  are  the  following: 


fo(0)  "  0 
fo(0)  -  0 
fo(0)  “  0 
fo(l)  “  1 
fo(i)  “  0 
fo(l)  -  0 


fo(l)  -  0 


fl(0)  -  0 
fl(0)  -  -1 
fl(0)  -  0 
fi(i)  -  0 
fl(l)  -  0 
f'ld)  -  0 


f'ld)  -  0 


fgCO)  -  0 

f2(0)  -  0 

f2(0)  -  -1 

fod)  -  0  (12) 

f2d)  ”  0 

f2(l)  “■  0  (used  only 

for  fifth- 
and  sixth- 
degree 

...  polynomials) 

f2(l)  “  0  (used  only 
for  sixth- 
degree 

polynomials) 


where  the  dot  denotes  dif ferentiation  with  respect  to 
Physically,  these  conditions  correspond  to  the  zero  slip  con¬ 
dition  at  the  wall,  satisfy  the  monentum  equation  at  the  wall 
through  the  second  and  third  derivatives  of  the  velocity  profiles, 
and  insure  a  smooth  transition  to  the  external  inviscid  flow.  As 
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a  result  of  the  use  of  these  boundary  conditions, 
become : 


f 

o 


fl 


f 


2 


r  1/3  (4ti  -  or 

<  1/3  (5Tj  -  +  3ri^)  or 

2r|  -  5ri  +  Gt)  -  2r\ 

1/6  (2^  -  or 

-<  1/4  (ri  -  2ri^  +  2^"^  -  or 

1/10  (2Ti  -  -f-  lOrj'^  -  lOri^  +  3^®) 

f  1/18  (ri  -  3^^  +  2Tj^)  or 
j  1/36  (rj  -  eff^  +  -  3t)^)  or 

1/60  Cn  -  10^^  +  2oii‘^  -  i5n^  +  4^®) 


the  polynomials 


(13) 


9.  Total  enthalpy  profiles  are  used  primarily  herein  whereas 
in  Part  I  (reference  (a))  static  enthalpy  profiles  were  used. 
While  the  choice  of  total  enthalpy  profiles  results  in  some 
simplification  of  the  mathematical  expressions  involved  it 
does  impose  a  limitation  on  the  method  which  is  discussed  in 
Paragraphs  20  through  22.  The  form  of  the  total  enthalpy  pro¬ 
file  is  the  same  as  that  of  the  static  enthalpy  profile;  namely; 


,  //  ,  , 


//o 


H 


o 


H2> 


are  defined  by  equations 


where  and  \ 

17  and  18.  The  physical  boundary  conditions  are  at  y  “0, 

H  “  Hw,  and  the  second  and  third  derivatives  of  the  enthalpy 
profiles  satisfy  the  energy  equation,  and  at  y  -  oo,  H  - 
and  the  first  three  derivatives  of  the  enthalpy  profiles  are 


zero . 
and , 

The 
hence , 

corresponding  transformed  plane  boundary 
the  polynomials,  are  therefore: 

conditions 

Go(0) 

“  1 

Gi(0) 

- 

0, 

G2  (0) 

n 

0 

Go(0) 

=  0 

Gi(0) 

rs 

-1 

G2(0) 
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0 
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-  0 
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- 
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G2(0) 
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1 
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m  0 

G^d) 
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God) 
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0 
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- 
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“  0 
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^  ,  _  -4  _5  _6 

1  “  +  2t]jj 

G],  (1h)  -  1/10  (s'hjj  -  +  lonj  -  loti^  +  3'^i®)  (le) 

^2  "  1/60  (-Tljj  +  10%  “  20%  +  15t|j|  - 


10.  The  equations  (11),  (12),  (14),  and  (15)  result  in  the 
following  definitions  of  the  profile  parameters 


The  expressions  for  the  profile  parameters,  or  form  factors, 
are  then  derived  for  the  zero  pressure  gradient  case  by  direct 
differentiation  of  the  transformed  momentum  and  energy  e- 
quations  (equations  (5)  and  (6)).  These  expressions  are: 

G  /^y.ojA/oY, 


dYx.oJ  A/ cf  fy^ 
/\*  J  L 


Gyy,o/- 


y-Pr^)[ i/A  *fyr,  o)  7  2 
//<.  ■ 


(18a) 


fyj^oy:oj[iyA*ryy,o)j^  [fP^^ Pr^  yj 
y  <S  //Vz 

^  Pw-^  w  f/l 


where  the  Y's,  which  are  determined  from  real-gas  property 
values  and  hence  are  functions  of  pressure  and  enthalpy,  are 
defined  as: 
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dh 

(18b) 

3  /‘2 

on 

Y_ 


and  A’*'  is  the  ratio  of  the  thermal  and  viscous  thicknesses, 
i  .3.  ,  A*  -  A. 


THERMODYNAMIC  AND  TRANSPORT  PROPERTIES 

li.  Since  the  real-gas  characteristics  of  the  present  method 
are  dependent  on  the  air  property  values,  only  the  most  recent 
data  on  the  properties  of  air  in  equilibrium  dissociation  have 
been  used.  Those  values  have  been  generated  by  the  Thermo¬ 
dynamics  Section  of  the  National  Bureau  of  Standards  and  were 
obtained  from  references  (g),  (h),  and  (i).  The  density,  as 
used  in  the  computations,  was  taken  inversely  proportional  to 
the  temperature  at  constant  pressure  for  temperatures  up  to 
2520OU.  For  temperatures  greater  than  2520OR,  the  density 
was  then  determined  as  a  function  of  pressure  and  enthalpy  by 
curve  fits  to  the  data  of  reference  (h)  as  developed  by  .Schmid, 
reference  (j).  In  addition  to  an  accurate  knowledge  of  the 
densi ty  it  was  necessary  to  have  reliable  values  of  the 
viscosity  and  Prandtl  number  so  that  the  derivatives  of  the 
quantities  Inpp  and  Inpp/Pr  could  be  formed.  The  values 
given  in  reference  (i)  were  considered  adequate  for  tempera¬ 
tures  below  2700*^;  however,  for  temperatures  between  2700OR 
and  6000OR,  new  constant-pressure  transport-property  tables 
were  generated  by  Klein,  reference  (g) ,  presenting  the  data 
at  small  temperature  increments  to  facilitate  differentiation. 
Since  the  derivatives  of  the  property  values  are  needed  only 
at  wall  temperatures,  the  tables  of  reference  (g)  were  termi¬ 
nated  at  6000°K  which  is  considered  more  than  adequate  for 
practical  wall  temperatures.  The  only  transport  property 
information  needed  at  higher  temperatures  was  the  viscosity 
for  the  computation  of  the  stream  Reynolds  numbers.  Viscosities 
at  temperatures  above  GOOQOR  were  determined  by  the  interpolation 
scheme  of  reference  (f) . 
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12.  The  need  for  accurate  determination  of  the  property 
derivatives  (i.e.,  equation  (18b))  is  apparent  from  their 
contribution  to  the  velocity  and  enthalpy  form  factorsj  equation 
(18a).  In  fact,  under  some  conditions,  inaccuracies  in  these 
derivatives  could  result  in  the  loss  of  solutions  as  will  be 
discussed  in  the  section  on  limitations.  In  an  effort  to 
obtain  maximum  accuracy  of  the  derivatives,  the  data  of  refer¬ 
ences  (g)  and  (i)  were  smoothed  by  the  numerical  procedtires 
described  in  reference  (k) ,  numerically  differentiated  by  the 
methods  outlined  in  Chapter  IV  of  reference  (1)  and  re-smoothed 
to  obtain  the  first  derivative.  This  procedure  was  repeated, 
using  the  results  of  the  first  differentiation,  to  obtain 

the  second  derivative.  It  was  then  necessary  to  make  slight 
additional  adjustments  to  assure  the  proper  agreement  between 
the  derivatives  and  the  original  quantities.  The  values  used 
are  presented  in  Figure  1, 

SOLUTION  OF  THE  SYSTEM  OF  EQUATIONS 

13.  The  basic  system  of  equations  to  be  solved  consists  of 
equations  (7),  (8),  (9a),  (11),  (14),  and  (18a),  Examination 
of  the  equations  indicates  that  a  solution  may  readily  be 
obtained  if  the  profile  parameters  are  assumed  to  be  inde¬ 
pendent  of  X  (also  Hq,  h^,  U,  and  P  are  constant).  This  is 
physically  equivalent  to  the  assumption  of  similar  profiles 
which  is  consistent  with  other  zero  pressure  gradient  analyses 
(e.g.,  references  (c) ,  (d) ,  (e) ,  and  (f)). 

14.  As  a  result  of  the  above  assumptions,  equations  (7)  and 
(8)  may  be  integrated.  Using  the  definitions  of  shear  stress 
and  heat  transfer  as : 

e,s 


where  the  subscript  w  refers  to  quantities  at  the  wall,  the 
integration  gives,  (for  0-E^Oatx-O): 

^2_  ZfP^fJw  _£  (21) 


(19) 


(20) 
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(22) 


15.  We  now  develop  equations  (9a)  so  that  the  ratios  6/6 
and  E/il  may  be  expressed  in  terms  of  the  profile  parameters 
and  the  bound ary- layer  thickness  ratio  A*.  This  is  ac¬ 
complished  by  inserting  the  polynomials  for  the  velocity  and 
total  enthalpy  in  equations  (9a)  and  performing  the  indicated 
integration.  The  results  are  expressed  in  terms  oi's  and 
Hi  (A*)  where  the  a^^’s  are  constants  the  Hi{A*)'s  are 
polynomials  in  A*,.  The  constants  and  the  coefficients  for 
the  polynomials  are  presented  in  Tah)le  I  for  fourth-,  fifth-, 
and  sixth-degree  velocity  profiles  and  for  A*  greater  or  less 
than  one.  The  required  ratios  are: 


-  % 


-  2^3  A,  A; 


(23) 


_  -  .  5'  ■A'2 

-f-  //  /<2i^  A,  -h//^ •/  A'y/hV  X 

The  additional  quantity: 


(/ 


=  /-  C*/o  X. 


:24) 


(25) 


is  developed  from  equation  (9b)  in  a  similar  manner;  however, 
the  values  of  ^/A  must  be  determined  by  direct  integration 
after  the  values  of  and  are  known. 


‘Hi 


H2 
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16e_  If  ejguations  (11)  and  (14)  are  differentiated  with  respect 
to  T)  and  T)jj  and  given  in  terms  of  the  wall  values  (i.e,,  rj  « 

■«  0)  we  have 

•  •  *  • 

(26) 


4  ^.fo)  *  ;{^G^  ^oj  (27) 

These  two  equations,  equations  (21)  through  (24),  and  the 
equation  (18)  now  independent  of  x,  constitute  the  systems  of 
equations  which  were  actually  solved. 

17.  The  method  of  solution  was  iterative  and  consisted  of  the 

following  steps.  First  values  of  X,  and  were  assumed, 

X  ^  •  2 

was  then  computed  from  equation  (26)  and  next  G^/A*,  X^  /A** 

and  Xjj  /A*^  were  determined  from  the  first,  third,  and  iourth 

equation  (18) .  These  quantities  were  inserted  in  equation  (27) 
and  a  cubic  in  A*  was  formed.  The  root  of  the  cubic  to  be  used 
was  decided  on  the  basis  of  physical  reasoning  and  then 
Xjj^,  and  Xfjg  were  computed.  Now  all  quantities  of  interest 

were  available  and  it  was  necessary  to  check  the  assumed  values 
of  Xj^  and  X2.  The  second  equation  (18)  gave  a  direct  check 
on  X2;  however,  it  was  necessary  to  combine  equations  (21) 
through  (24)  to  check  X^  in  terms  of  the  computed  When 

the  differences  between  the  computed  quantities  and  the 
quail  lilies  as  determined  by  the  checking  equationfs  was  less 
than  0.001  percent,  '.ne  .solution  was  considered  satisfactory. 

If  the  checking  equation  was  not  satisfied,  X^  and  X2  were 
incremented  by  both  positive  and  negative  quantities  and  the 
iteration  process  was  repeated  until  satisfactory  results 
were  achieved. 


18.  After  the  system  of  equations  was  solved,  the  values  of 
d*  and,yi.were  obtained  from  equation  (25)  and  equation  (9b), 
respectively,  l^ext,  the  inverse  transformations,  equations 
(10) ,  were  applied  to  obtain  the  boundary-layer  thicknesses 
in  the  physical  plane  as  a  function  of  the  square  root  of  x; 


I’he  shear  stress  and  heat 
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transfer  at  the  wall  may  now  be  expressed  as  function  of  1/  x 
by  equations  (19)  and  (20).  If  the  local  skin-friction  coef¬ 
ficient,  Cf,  is  desired  it  may  be  derived  from  equation  (21) 
as : 


^  /f 


z 


9/^ 


(28) 


where  O/n  is  given  by  equation  (23)  and  Cf  is  the  mean  friction 
coefficient.  Similarly,  if  a  Stanton  number, 

St  -  _  - 

is  desired,  equation  (29)  may  be  developed  Lo  give: 


(29) 


where  E/A  is  given  by  equation  (24)  and  hj^  is  the  recovery 
enthalpy . 

19,  When  computing  the  recovery  enthalpy  as  defined  by  the 
relations : 


and 

fi.-h,  . 

some  investigators,  (e.g.,  references  (d)  and  (e))  have 
found  that  R  is  a  function  of  Mach  number  and  wall--to-3tream 
temperature  ratio.  The  heat- transfer  results  presented  herein 
were  calculated  using  total  enthalpy  profiles  and  are  based 
‘on  a  unity  recovery  factor  to  avoid  ambiguity  as  to  the  choice 
of  recovery  factor  used  in  the  computations.  The  heat- transfer 
results  of  reference  (f)  have  also  been  adjusted  to  a  unity 
recovery  factor.  For  purposes  of  comparison  some  recovery 
factors  have  been  computed  for  the  case  of  an  insulated  wall 
using  only  static  enthalpy  profiles.  Difficulties,  discussed 
in  Paragraph  21,  were  encountered  when  the  total  enthalpy 
profiles  were  used  to  derive  the  insulated  wall  recovery 


(30) 


(31) 
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factor.  The  recovery  factor  for  the  static  enthalpy  profiles 
is : 

R  ^  O.Z  (a.'' (32) 


where  f^  is  found  from  equation  (26)  with  equal  to  zero. 
The  quantity  A*  is  found  from  a  form  of  equatiftn  (24)  derived 
with  a  static  enthalpy  profile  and  with  E/A,  and  Xjj 

equal  to  zero.  ^ 


LIMITATIONS  OF  METHOD 

20,  The  one-parameter  velocity  profile  representation,  as 
described  in  Part  I,  reference  (a),  results  in  a  limitation 
to  the  range  of  applicability  of  the  present  method.  This 
limitation  occurs  because  under  most  physically  realistic 
flow  conditions  the  absolute  value  of  the  profile  parameter 
cannot  exceed  12.  Physically,  the  existence  of  profile  param¬ 
eters  less  than  minus  12  indicates  that  flow  separation  would 
occur  and  if  the  parameter  exceeded  plus  12  it  indicates  that 
velocity  overshoot  occurs  within  the  boundary  layer.  While 
the  latter  case  is  possible  under  some  conditions  of  heat 
transfer  it  is  normally  not  probable  for  cases  with  heat 
transfer  to  the  surface  and,  therefore,  is  a  limitation  on 
physical  grounds.  Actually  the  significant  root  of  the  equation 
for  the  velocity  profile  parameter  became  imaginary  in  the  one 
parameter  method  after  its  value  reached  12.  The  lira^itations 

of  flow  separation  and  velocity  overshoot  are  general  and  also 
apply  to  the  two-parameter  velocity  profile  representation. 

In  this  case,  the  limiting  values  of  X]_  and  X2  are  found  by 
equating  the  derivative  of  the  velocity  profile  to  zeio,  l.e., 
f(x,Ti)  -  0,  which  gives  a  linear  relation  between  the  variables 
X^  and  X2  for  a  fixed  value  of  n-  Then  q  is  eliminated 

from  the  family  of  such  linear  relations  yielding  an  envelope 
whose  interior  represents  the  allowable  values  of  X^^  and  X2. 

The  envelopes  vary  with  the  degree  of  the  polynomial  and  are 
presented  in  Figure  2  for  fourth-,  fifth-,  and  sixth-degree 
polynomials.  It  is  noted  that  the  higher  degree  polynomial 
yields  a  larger  envelope.  While  these  envelopes  were  not 
exceeded  for  any  of  the  flat-plate  cases  presented,  they  are 
useful  as  an  indication  of  the  validity  of  solutions  obtained 
in  more  general  cases  such  as  with  pressure  gradient. 

21.  In  a  somewhat  similar  manner,  the  choice  of  sixth-degree 
total-enthalpy  profiles  imposes  a  limitation  on  the  values  of 
the  wall  heat  transfer  which  may  be  obtained.  This  limitation 
results  in  an  inability  to  represent  the  zero  heat-transfer 
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total-enthalpy  profile  and  implies  that  calculations  near 
zero  heat  transfer  may  be  suspect.  In  this  case,  A 

zero  and  G  (^„)  •*  0  yields  a  condition  when  •»  0  such 
n  ,  11 

that  no  further  values  of  G  (^pj)  “  0  are  possible  for  iijj 

between  zero  and  one^  An  additional  zero  derivative,  however, 
is  necessary  in  the  range  to  represent  properly  the  zero 
heat-transfer  total-enthalpy  profile.  Fortunately,  this 
limitation  does  not  exist  when  either  the  pressure  gradient 
or  wall  enthalpy  gradient  is  non-zero. 

22.  Another  limit  which  is  applicable  to  the  zero  pressure 
gradient  case  is  brought  about  by  unfavorable  combinations 
of  property  values  (i.e.,  equation  (18b)),  and  flow  conditions. 
By  combining  the  first  two  eqiiations  of  (18a)  with  equation 
(26),  a  quadratic  equation  may  be  formed  in  terms  of  A  , 

When  the  resulting  equation  is  solved  for  1  ^  the  discrimi¬ 
nant  may  become  negative  and,  hence,  lead  to  complex  values 
of  '  This  limit,  which  varies  with  the  degree  of  the 
velocity  polynomial,  is  not  violated  if  the  following 
inequalities  are  satisfied. 

Fourth  degree  ~/^/^  ( 2/i  -  2.  ^ ^  O  (33) 


n  /O  r  ^  S 
Fifth  degree  -p/j  r  — ~ 
/  /  Z  ‘w  4 


■2>)^<POf3 


■o 


(34) 


Sixth  degree 


(35) 


where  = 


>r  yp 


(36) 


Since  and  B2/ (Hq  -  b^)  are  functions  of  wall  temperature 
and  pressure  it  is  difficult  to  generalize  with  respect  to 
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the  above  inequalities.  Figure  3  has  therefore  been  plotted 
for  0.1  atmosphere  pressure  and  presents  the  value  of  the 
left-hand  member  of  the  inequality  (designated  as  G  (Si,S^) 

as  a  function  of  h^/RTQ,  where  RTq  -  33.7098  BTU/lb,  for 
constant  values  of  f^  and  -  hj^.  Since  must  be  positive 

and  generally  lies  between  one  and  two,  this  presentation  gives 
a  qualitative  guide  to  the  combination  of  conditions  which 
should  fail  to  give  solutions.  The  curves,  which  are  ex¬ 
ceedingly  irregulai*  because  of  the  real-gas  properties  used, 
i.e,.  Figures  Ic  and  Id,  indicate  that  the  criterion  should  be 
examined  on  a  per  case  basis.  It  was  also  noted  that  changes 
in  Y  which  were  smaller  than  the  expected  accuracy  of  the 

numerically  determined  second  derivatives  could  change  a  case 
from  violating  to  satisfying  the  criterion.  From  this  it 
must  be  concluded  that  still  further  accuracy  is  needed  in 
determining  the  property  derivatives,  i.e.,  the  equation  (18a) 
quantities.  Figure  3  shows  in  addition  that  excessively 
large  values  of  (Hq  -  hi)  and  small  values  of  ±2  have  a  greater 

probability  of  failing  to  yield  solutions.  One  may  also 
observe  that  for  a  perfect  gas  with  constant  Prandtl  number 
and  constant  coefficient  of  specific  heat  and  viscosity  described 
by  the  Sutherland's  law  the  value  of  fi  is  always  positive  and 

hence,  the  criterion  is  satisfied  in  almost  all  practical  cases. 

RESULTS  AND  DISCUSSION 

23.  It  was  pointed  out  in  the  section  on  limitations*  that 
the  use  of  total-enthalpy  profiles  precluded  investigation 
of  the  zero  heat  transfer  (i.e.,  insulated  wall  case).  In 
order  to  obtain  a  skin  friction  Mach  number  relatioi-.sbln  for 
zero  boat  transfer,  a  limited  number  of  computations  were  made 
with  static  enthalpy  profiles  and  sixth-degree  velocity  pro¬ 
files.  The  results,  in  terms  of  Cf/l'vs.  Mx,  are  pres¬ 
ented  in  Figure  4.  In  Figure  5  the  corresponding  insulated 
wall  enthalpy  recovery  factor  is  presented.  These  compu¬ 
tations  were  made  using  equation  (32)  and  also  by  using  the 
square  root  of  the  Prandtl  number  taken  at  the  derived  wall 
temperature  in  zero  heat  transfer.  The  present  values  are 
compared  with  those  of  references  (d)  and  (e) .  From  Figure 
4a  it  is  seen  that  the  agreement  in  Cf  /  R  is  reasonable  at 


♦  Because  of  the  indicated  limitations  all  results  are  given 
in  the  figures  with  symbols  used  to  represent  actual  computed 
values.  In  some  cases  the  curves  have  been  faired  through 
regions  where  solutions  are  not  attainable. 
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all  stream  Mach  numbers  with  slight  discrepancies  occurring 
at  Mach  number  near  one.  Correspondingly  good  agreement  is 
noted  in  Figure  5  for  the  insulated  wall  enthalpy  recovery 
factor  where  the  greatest  discrepancy  between  the  methods 
is  in  the  order  of  five  percent.  As  has  been  noted  by  other 
investigators,  the  enthalpy  recovery  factor  curve  varies 
closely  as  the  square  root  of  the  wall  Prandtl  number.  These 
results  indicate  that  if  investigations  are  to  be  conducted 
by  the  two  parameter  integral  method,  use  of  static  enthalpy 
profiles  should  give  reasonable  results  for  cases  near  zero 
heat  transfer,  whereas,  the  use  of  total-enthalpy  profiles 
would  fail  to  give  satisfactory  values. 

24.  Since  the  case  of  equal  wall  and  stream  temperatures 
reduces  to  the  zero  heat-transfer  case  at  zero  Mach  number 
it  was  considered  a  desirable  case  for  comparing  the  results 
obtained  by  using  velocity  profiles  of  varying  detgree.  Fourth-, 
fifth-,  and  sixth-degree  velocity  polynomials  were  used  as 
well  as  a  sixth-degree  total  enthalpy  presentation.  The 
computed  values  of  Cf  \/  Re  are  shown  in  Figure  4  as  a  function 
of  Mach  number.  As  might  be  expected  the  higher  degree  pro¬ 
file  polynomials  give  the  better  agreement  with  the  less 
approximate  methods.  The  sixth-degree  presentation  is  within 
less  than  one  percent  of  the  Blasius  value  at  zero  Mach  number. 
As  the  Mach  number  is  increased,  however,  the  results  of  the 
methods  (i.e,,  references  (c)  and  (e) ,  and  the  present  method) 
diverge  considerably,  differing  by  as  much  as  20  percent  at 
Mach  number  20.  As  will  be  seen  later,  this  difference  is  not 
typical  of  most  cases  but  may  be  interpreted  as  an  Indication 

of  a  difficlency  in  the  present  method  when  computing  cases 
where  the  wall-to-stream  static  enthalpy  ratio  is  in  the  order 
of  one  and  the  Mach  number  is  large. 

25.  The  methods  of  references  (c),  (d),  and  (e)  do  not  include 
the  effects  of  air  dissociation  on  the  thermodynamic  and  trans¬ 
port  properties  and,  hence,  comparisons  with  the  two- thickness 
integral  method  in  regions  where  dissociation  is  present  are 

at  best  qualitative.  The  recent  work  of  Wilson,  reference 
(f),  does  present  an  opportunity  for  quantitative  comparison 
in  the  dissociation  region.  Following  Wilson  the  skin-  friction 
is  presented  in  Figures  6  and  7  as  a  function  of  wall-to-stream 
static  enthalpy  ratio  for  constant  values  of  the  parameter 

U^/2h^  at  stream  enthalpies  of  93.24  and  466.3  BTU/lb  (i.e., 

3000  and  15,000  BTU/siug),  Both  fourth-  and  sixth-degree 
velocity  polynomials  and  total-enthalpy  profiles  were  used  to 
obtain  the  results  which  are  presented  and  again  the  agreement 
is  considerably  more  favorable  when  the  sixth-degree  polynomials 
are  used.  The  results  obtained  with  the  sixth-degree  velocity 
profiles  are  generally  within  approximately  five  pereeat  of 
the  reference  (f)  values  with  the  notable  exception  of  the 
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region  of  low  wall-to-streaia  enthalpy  ratio.  The  reason  for 
this  difference  Is  not  clear  nt  present  but  is  believed  to 
be  associated  with  the  occurrence  of  an  inflection  point  in 
the  transformed  plane  velocity  profiles  which  is  uniquely 
fixed  in  q  by  the  polynomial  approximation.  The  transformed 
plane  inflection  point  itself  appears  to  be  necessary  to 
accommodate  the  large  static  enthalpy  gradients  which  are 
present  in  these  highly  cooled  wall  cases, 

26.  In  Figure  8  the  present  results  ai‘e  compared  with  those 
of  reference  (f)  on  the  basis  of  a  modified  Reynolds  analogy 
parameter,  ^ 

13.L3r^  j 


The  Stanton  number  for  the  present  results  were  determined 
using  a  unity  recovery  factor.  It  therefore  necessary 

to  alter  tlie  reference  (f)  results  to  include  a  unity  recovery 
factor  so  that  a  common  basis  of  comparison  was  achieved. 

Here  again,  the  sixth-degree  velocity  profiles  appear  to  give 
a  better  representation  of  the  physical  picture.  With  respect 
to  heat  transfer,  however.  Figure  9  shows  good  agreement 
between  the  two  methods  which  implies  that  the  differences 
indicated  in  the  modified  Reynolds  analogy  factor  are  due 
to  the  previously  discussed  differences  in  the  skin  friction. 
With  respect  to  the  two-thickness  integral  method,  it  appears 
that  the  use  of  either  the  fourth-  or  sixth-degree  velocity 
profile  representation  give  essentially  identical  values  of 
the  heat  transfer. 

27.  A  limited  investigation  of  the  effect  of  pressure  level 
on  the  skin-friction  coefficients  and  modified  Reynolds 
analogy  factor  was  conducted.  The  results  presented  in 
Figures  10  and  11  do  not  show  any  striking  features  except 
possibly  a  larger  effect  of  wall-to-stream  enthalpy  ratio 

at  the  larger  pressure  levels. 

28,  As  a  final  comparison  between  the  results  of  the  reference 
(f ) ,  and  the  present  method,  velocity  and  static  enthalpy 
profiles  in  the  physical  plane  are  presented  in  Figure  12. 

The  case  presented  is  for  a  wall-to— stream  enthalpy  ratio  of 
one  and  as  previously  discussed  (Paragraph  25)  represents  the 
area  of  greatest  divergence  between  the  two  methods.  It  is 
interesting  to  note  that,  in  the  lower  half  of  the  boundary 
layer,  the  profiles  are  quite  similar.  While  it  is  difficult 
to  generalize,  this  may  indicate  that  quantities  such  as  heat 
transfer  and  skin  friction  which  are  dependent  on  derivatives 
of  the  profiles  at  the  wall  could  be  reasonably  well  repre¬ 
sented. 
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CONCLUSIONS 

29.  In  general,  it  may  be  concluded  that  the  two-thickness 
integral  method  utilizing  a  sixth-degree  polynomial  velocity 
profile  representation  gives  reasonable  results  when  compared 
with  the  "exact"  methods. 

30.  The  present  investigation  has  also  revealed  limitations 
of  the  two-parameter  polynomial-approximation  procedure  most 
of  which  must  be  examined  on  a  per  case  basis.  Among  the 
limitations  are  the  following: 

a.  The  possibility  of  loss  of  solutions  due  to  un¬ 
favorable  combinations  of  thermodynamic  and  transport  property 
derivatives . 

b.  The  failure  of  total  enthalpy  profile  polynomials  in 
cases  of  zero,  or  near  zero,  heat  transfer. 
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Officer,  Office  of 
Naval  Res.  Br.  Office 
346  Broadway 
New  York  13,  N.  Y. 

Sandia  Corporation 
Sandia  Base 

Albuquerque,  New  Mexico 
1  Attn:  Mr,  Alan  Pope 

1  Defense  Metals  Inf.  Center 
Battelle  Memorial  Institute 
505  King  Avenue 
Columbus  1,  Ohio 

Commanding  General 
Army  Rocket  and  Guided 
Missile  Agency 
Redstone  Arsenal,  Alabama 
1  Attn:  John  Morrow 

National  Bureau  of  Standards 
Washington  25,  D.  C. 

1  Attn;  Dr.  G.  B,  Schubauer 

Cornell  Aero.  Laboratory 
4455  Genesee  Street 
Buffalo,  N.  Y, 

1  Attn;  Dr,  Gordon  Hall 

General  Applied  Sciences 
Laboratories,  Inc, 

Merrick  and  Stewart  Avenues 
East  Meadow,  New  York 
1  Attn:  Mr.  Robert  Byrne 

Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena  3,  California 

1  Attn:  I.  R.  Kowlan,  Chief 

Reports  Group 

2  Attn;  Dr.  L.  Jaffee 
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EXTERNAL  DISTRIBUTION  LIST  (SP) 

No.  of 
Copies 

Los  Alamos  Scientific  Lab. 

P.-  O.  Box  1663 
Los  Alamos,  New  Mexico 
1  Attn:  Dr.  D.  F.  MacMillan 

N-1  Group  Leader 

Institute  for  Defense  Analyses 
Advanced  Research  Projects  Agency 
Washington  25,  D.  C. 

1  Attn:  Mr.  W.  G.  May 

General  Sciences  Branch 

Kaman  Aircraft  Corporation 
Nuclear  Division 
Colorado  Springs,  Colorado 
1  Attn:  Dr.  A.  P.  Bridges 

U.  S.  Atomic  Energy  Commission 
P.  0.  Box  62 
Oak  Ridge,  Tenn. 

1  Attn:  TRI : NLP: ATD: 10-7 

1  Sandia  Corporation 
Livermore  Laboratory 
P,  0,  Box  969 
Livermore,  California 

United  Aircraft  Corporation 
Research  Laboratories 
East  Hartford  8,  Conn. 

1  Attn:  Mr.  II,  J.  Charette 

1  Attn;  Mr.  H.  Taylor 
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